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SUMMARY 

A one-dimensional analysis of nonconstant area and nonconstant pressure burning 
and mixing processes is carried out using a prescribed pressure-area relation. The duct 
area ratio is advanced as an important parameter that affects the existence of flow solu- 
tions and the loss or buildup of the stagnation pressure ratio of the process. The critical 
condition of thermal choking associated with Mach number of unity for constant area duct 
is shown to assume other Mach values for various nonconstant area ducts. The concepts 
evolved are applied to idealized supersonic combustion ramjets and rocket-ramjet hybrid 
engines. 

The nonconstant area combustion in a supersonic combustion ramjet results in a 
better specific impulse. Also, the separate mixing and burning, each with its appropriate 
duct configuration, leads to better performance than simultaneous mixing and burning in 
the case of the augmented rocket. 


INTRODUCTION 

The analysis of combustion or mixing processes in jet engines is usually confined to 
the cases of constant cross-sectional area or of constant gas pressure. References 1 to 
3 are representative of the scope of the work in this field performed under the previously 
mentioned constraints. The present analysis is an attempt to treat the cases going be- 
yond these confining stipulations. The burners and mixers thus evolved are applied to 
supersonic combustion ramjets and rocket- ramjet hybrid engines. In order to preserve 
the overall view of the processes involved, a one-dimensional treatment is carried out 
that assumes ideal gases with specific heats that do not change with temperature, no 
frictional penalties, and hydrogen stoichiometric burning in the primary or secondary 
flow. Complete mixing is assumed. 

First, the supersonic burning is dealt with under the conditions of a variable geome- 



try duct, and the requirements for its optimum performance are outlined. A critical con- 
dition, analogous to thermal choking at Mach 1 for the constant area duct, has to be re- 
cognized: this condition occurs at a Mach number that is not necessarily unity. The 
burning which takes place with its Mach number decreasing toward the critical is refer- 
red to as supercritical burning, while the burning which proceeds with its Mach number 
increasing toward the critical is referred to as subcritical burning. Thus, the concept of 
supersonic combustion may be replaced by the concept of supercritical combustion. 

There are also duct configurations where the Mach number hardly changes due to burning, 
however extensive this burning is. 

Next the mixing is investigated by exploring (1) conditions for optimum buildup of the 
stagnation pressure of the secondary flow and (2) the option of a separate mixing and burn- 
ing as compared to the simultaneous mixing with burning. The analyses are applied in 
the consideration of the performance of the ramjet and rocket-ramjet hybrid engines, 
and the ways of arriving at optimum duct configurations are illustrated. The appendix 
contains the derivation of pertinent equations. 


SYMBOLS 



F 

G 

G* 


e 

h 

J s 

J 

m 


cross-sectional area of duct 

area ratio of secondary inlet over primary inlet, A^/A^ 


duct area ratio, exit over inlet 

net thr ust 

(dynamic pressure)(A^ + A^) 


thrust coefficient, 


function of Mach number, M-j/ 1 +. ? - — 




function of Mach number and power index, e + yM 

.« _ o € (r- 1 >A' ( £ + yM 2 ) C * y ~ 1>/y 


function of Mach number and power index, G = 


X 


1 + 


hIm 2 


n 

dimensional constant, 32. 2 (lb mass - ft)/(lb - sec ) 
enthalpy, Btu/lb of gas 
specific impulse, sec 

mechanical equivalent of heat, 778 ft - lb mass/Btu 
mass flow per second, lb mass/sec 
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rrtj entrainment ratio, m^/mu 

M Mach number 

M*/l + M 2 

N function of Mach number and power index, — — 

e + yM^ 

o 

P pressure, lb/ft 

q heat added, Btu/lb of air 

R gas constant, ft-lb/(lb mass)(°R) 

S entropy, Btu/(lb)(°R) 

T absolute temperature, °R 

X parameter, 1 + 2 i-l_i m^ 

2 

V velocity, ft/sec 

a parameter, sin a = N/N. 

y ratio of specific heats, C /C 

r “ 

A ratio of stagnation enthalpy of secondary flow h° over stagnation enthalpy of pri- 
mary flow h? 

e power index in L. Crocco power area - pressure relation, P2/ p l = 
p density, (lb mass) /ft 3 

Subscripts: 
av average 

b pertaining to burning 

c value of function at critical point 

j primary flow inlet 

m end of mixing 

0 referring to ambient conditions 

1 secondary flow inlet 

2 end of burning 

3 end of nozzle expansion 
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Superscript: 


o stagnation condition 


PRESSURE - AREA RELATION 


In solving the duct flow problems of heat addition or mixing, the one -dimensional 
treatment of the three conservative equations in their integrated form is found very use- 
ful for many flow regimes. A convenient assumption of the constancy of a certain quan- 
tity throughout the flow (e. g. , the pressure or duct cross-sectional area) is commonly 
made. This permits the evaluation of the wall force integral in the momentum equation 
and hence the retention of the simple integrated form. 

In order to enlarge the scope of the flow solutions, it would be desirable to have 
something else besides constant area or constant pressure cases. In general, there are 
many ways of expressing mathematically the variation of pressure with flow area. The 
one that has been chosen for this analysis because of its simplicity is due to L. Crocco 
(ref. 4) and is represented by 




e/(l-0 



Figure 1. - Pressure-area relation. 



The pressure at any axial station is related to flow 
area through the single parameter e. It is to be rec- 
ognized that this arbitrary assumption as to the vari- 
ation of the pressure inside the duct in no way condi- 
tions the length of the duct. Therefore, the area of 
the burner or mixer could be varied with length in such 
a manner that the existing finite rates of heat release 
or mixing yield the assumed variation of pressure 
with A. The convenience of this assumption becomes 
apparent not only by the ability to retain a simple inte- 
grated form but by the inclusion of constant pressure 
case with e = 0 and constant area case with e = 1. 
Thus, the introduction of the e relation leads to some 
generalization and extension of the two more common 
cases. 

Figure 1 presents the variation of the pressure 
ratio against the area ratio for several e’s. All pos- 
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sible variations are exhausted by positive and negative e' s. 

It is of interest to study this generalized process in more detail - at first in combus- 
tion as an extrapolation of the known Rayleigh process and later on in mixin g. 


COMBUSTION 


The process of heat addition to a flow in a duct of a constant cross-sectional area can 
be described in terms of the familiar Rayleigh curve in the h - S diagram shown below 
in sketch (a). 


The extreme entropy point on the curve is a critical point 
depicting a thermal choking condition. The upper branch of the 
curve is associated with a subsonic flow, which with heat addi- 
tion has its Mach number increasing, and the lower branch is 
associated with supersonic heat addition flow, with Mach num- 
ber decreasing. The Mach number at the critical point is unity. 
Heat addition in excess of that required to achieve the critical 
or choking terminal Mach number will force an adjustment of 
the flow. The initial Mach number will be changed to a magni- 
tude that is consistent with the amount of heat input. 

It is fairly straightforward to show that the least stagnation pressure loss associated 
with a given amount of supersonic burning occurs in a process that reaches the critical 
condition. 




Generalized Rayleigh Process 

With the introduction of a nonconstant- 
area combustion process of the type specified 
by the Crocco e relation, the enthalpy - 
entropy diagram assumes the form shown in 
sketch (b). 

When one starts from the same initial 
value of the enthalpy, different e's trace out 
different curvature lines. The upper branches 
of the lines delineate a subcritical process 
where Mach number increases (finally reach- 
ing the critical condition) because of heat ad- 
dition, and the lower branches correspond to 
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a supercritical process where the Mach number is decreasing. One of the curves is the 
familiar Rayleigh line with e = 1. Again, it can be shown that the minimum stagnation 
pressure loss of the supercritical burning for fixed heat addition and e occurs when the 
entrance is chosen such that the critical condition is reached at the end. 


Critical Condition 

To obtain a complete description of the flow at the exit section of the burner when the 
initial Mach number Mj, the initial stagnation enthalpy hj, and the amount of heat added 
q are given, an energy balance equation is required which with the help of the other two 
conservation equations takes the form 


/^\ 2 =1+ q_ (7) 

W 

where 

M-|4~7-21^-i M 2 

N = — I 2 (6) 

e + yM 2 

(The equations are described in the appendix. The equation numbers remain the same. ) 
Equation (7) may be solved directly for Mg. When the function N is plotted against the 
Mach number M, figures 2(a) and (b) are obtained. In figure 2(a) for positive e’s, the 
highest points on the N curves define the critical conditions and result from the unique 
solutions to the previous heat addition equation. The critical Mach number is given by 



( 8 ) 


In the case of constant-area heat addition with e = 1, the critical Mach number is 1. 

For e < 1, this critical Mach number is subsonic, and for e > 1, the critical Mach num- 
ber is supersonic. A positive e varying from 0 to y/(y - 1) exhausts all possible values 
for M c from 0 to °°. In general, for a given value of e, two solutions or no solutions 
result as shown in sketch (c) when equation (7) is solved for Mg. Moderate heat addition 
at a given Mj will yield subcritical and supercritical Mg corresponding to a new Ng. At 
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maximum heat addition, there is a single 
solution at a critical Mg. In the case of 
constant-area heat addition with e = 1, the 
one-solution situation corresponds to the 
thermally choked flow. Greater heat addi- 
tion at a fixed yields cases where no 
solution (and hence no steady flow) is pos- 
sible. In the case of e = -y, the N curve 
increases rapidly in the neighborhood of 
M = 1, never quite reaching this asymptote. 
The maximum value of Ng in this case is 
dictated by stoichiometric fuel addition 
which maximizes q and hence Ng/Nj (see 
eq. (7)). This condition corresponds to a 
practically constant Mach number heat ad- 
dition (see sketch (d)). The duct configura- 
tion is such that a flow at a Mach number 
approaching unity remains near unity after 
heat addition. From figure 2(b) it is seen 
that with a negative e there is in general 
an asymptotic value of Mach number equal 
to V-e/y that the flow cannot quite reach; 
however, this Mach number can be ap- 
proached as close as desired whether on 
the subcritical or the supercritical side. 

In sketch (e) the convergence or divergence 
of the duct is shown dependent on whether 
the flow is supercritical or subcritical and 
on the value of e. 

Stagnation Pressure Ratio 

In order to be able to judge the merit 
of one particular e (in other words, the 
advantage of one particular duct configur- 
ation for the given initial flow conditions) 
and one particular heat addition, the in- 
fluence of e on the stagnation pressure 
loss during burning has to be considered. 


8 


I III 


— M ii iii! i*"* 1 nn i in ini"» i 



Figure 3. - Stagnation pressure ratio due to burning that changes flow from diffusion to critical Mach number. 



Figure 4. - Stoichiometric hydrogen combustion stagnation pressure ratio as function of diffusion Mach number 
for flight Mach number of 20. 


Figure 3 shows the stagnation pressure ratio achieved during burning initiated at 
and ending at M c . The stagnation pressure ratio of the process ending at any other Mach 
number is equal to the ratio of the pressure ratios corresponding to the initial and f inal 

P 2 p ? / P c 

Mach numbers since — = = — / — . The relation between M. and M„ is found by ap- 

P 1 P l/ P 2 

plying the N expression of the preceding section for a particular inlet stagnation enthalpy 
(function of flight Mach number) and heat addition (function of fuel type and fuel-air ratio). 
Figure 4 shows the stagnation pressure ratio for several € values and supercritical burn- 
ing for the specific case of flight Mach number Mq = 20 and stoichiometric combustion 
of hydrogen. It is seen that given a series of e ducts with the supercritical flows at the 
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same initial Mach number Mj and subjected to the same heat addition, the larger the e, 
the smaller will be the stagnation pressure loss. From sketch (e), where the shapes of 
the ducts are noted, this would imply, for example, that a convergent duct is superior to 
a constant area duct for the same Mj. 

Given a series of e ducts with the supercritical flows of the same stagnation enthalpy 
at different initial Mach numbers M^, such that each duct reaches a critical condition after 
the same heat addition, the larger the positive e, the larger will be the stagnation pres- 
sure loss (circled points in fig. 4). This implies that a constant area duct is superior to 
a convergent duct if a critical condition is reached in both ducts. (The shapes of duct are 
noted from sketch (e) . ) Thus diffusing (without losses) the flow sufficiently for the subse- 
quent heat addition in the constant area duct to reach a critical condition will yield a 
better stagnation pressure ratio than direct use of a critically convergent duct. A diverg- 
ing, critical, constant Mach number duct (e = - 1. 4) is still better. Best of all is the con- 
stant pressure case (e = 0); however, heat addition in this case requires the flow to pass 
smoothly from supersonic to subsonic. This condition may not be achievable in practice. 

In general, the critical e for any given Mj, for a given heat addition, and for ini- 
tial stagnation enthalpy, that secures the attainment of the critical endpoint, can readily 
be identified by a direct calculation as derived in the appendix (see eq. (15)). 


Optimization of Inlet-Combustor Combination 

In the design of an inlet-burner combination, the choice of e for the burner duct and 
the diffusion Mach number M^ will reflect directly on the overall stagnation pressure 
loss. For every e it is seen that the lowest permissible Mj minimizes the stagnation 
pressure loss due to heat addition. However, diffusion to low values of M^ generally 
causes greater pressure losses in the inlet diffuser. Selection of the optimum e and 
Mj must hence consider the combined inlet and combustor. 

The performance of a supersonic combustion ramjet inlet is a function of the inlet 
type, the flight Mach number Mq, and the design diffusion Mach number Mj. For illus- 
trative purposes, figure 5 shows the pressure recovery of a representative supersonic 
combustion ramjet inlet. It is based on the evaluation of the stagnation pressure of the 
flow following three oblique shocks of equal static pressure rise. 

If one has an inlet whose performance is a function of the diffusion Mach number Mj 
and a burner with a specific e, where the stagnation pressure ratio achieved is a func- 
tion of Mj, it is an easy matter for a given flight condition and heat addition to search 
and find an Mj that will result in an overall minimum stagnation pressure loss. 

Generally, the critical condition for burning is not reached. Diffusing the inducted 
flow of air to a very low Mach number Mj incurs a heavy penalty in stagnation pressure 
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Figure 5. - Inlet pressure recovery of a representative supersonic combustion ramjet for 
three oblique shocks of equal strength. 


recovery; however, this makes it possible for the burning to take place near the criti- 
cal point with its favorable stagnation pressure ratio. On the other hand, hardly diffus- 
ing the oncoming flow, and thus allowing a high pressure recovery, compromises the 
stagnation pressure ratio of combustion that of necessity occurs far off the critical point. 
Somewhere in between those two extreme conditions an optimum is usually found. 

If one uses the inlet data of figure 5 and the combustor data of figure 4, from which 
the stagnation pressure variation of hydrogen stoichiometric burning against Mj can be 
found, the effect of Mj on the overall pressure ratio can be shown as in figure 6 for 
Mq = 20 and a range of e. 

It is seen that the maximum stagnation pressure is reached for e = 2. 0 and 
= 8.2 under the condition of critical burning. Figure 7 shows the variations of the 
postcombustion Mach number Mg against the diffusion Mach number Mj. The postcom- 
bustion Mg hardly changes at large diffusion Mach numbers. 

Repeating the same procedure for a series of flight Mach numbers gives figure 8, 
where the optimum diffusion Mach numbers are noted. The signficance of the maximiza- 
tion of the stagnation pressure ratios product for the evaluation of the optimum specific 
impulse is seen from equations (12) and (13) in the appendix. The specific impulse evalu- 
ated at the ascertained Mj’s and e’s, which change with M Q , is shown in figure 9 for 
the Mq range of 10 to 25. This is compared with the performance of the constant e 
engines. It is seen that modest, but significant, improvements are offered by the 
optimum-e case as compared to the constant-area case. It is noted that the improvement 
would be magnified had a less idealized expansion nozzle been assumed. 

The e of the optimum variable e ramjet is seen to vary from e = 0. 5 at 
Mq = 10 to e > 2 at Mq = 25. Figure 10 describes in detail the e and the associated 
area ratio variation of the optimum ramjet for the assumed conditions. From a divergent 
burner duct at lower Mach numbers Mq, the specified duct shape changes to a constant 
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Figure 6. - Overall stagnation pressure ratio due to induction and combustion (hydrogen 
stoichiometric) for flight Mach number of 20. 



Figure 7. - Postcombustion Mach number as function of diffusion Mach 
number for flight Mach number of 20. 
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Figure 8. - Diffusion Mach number as function of flight Mach number. 



Flight Mach number, M 0 

Figure 9. - Ramjet specific impulse as function of flight Mach number. 
Hydrogen stoichiometric combustion; full expansion into ambient con- 
ditions. 
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area and later to a convergent one at 
higher flight Mach numbers. 

The indicated desirability of em- 
ploying a convergent combustor at high 
flight Mach numbers should be tem- 
pered by the following realization: as 
previously pointed out (fig. 4), best 
combustion performance is obtained if 
the Mach number entering the combus- 
tor is near unity and heat is then added 
with an € of 1 or less (corresponding 
to constant area or diverging area 
ducts), since this yields an exit Mach 
number near unity. However, diffusing to such a low Mj causes excessive inlet losses 
according to the assumed schedule of figure 5, so that the optimum Mj is considerably 
higher. In order to retain the benefit of a low exit Mach number, the calculations then 
indicate the use of a converging duct (provided by a high e) . The duct thus performs the 
function of a diffuser, but without the pressure losses associated with the assumed inlet 
(since the present calculations consider combustor stagnation pressure losses that arise 
from heat addition only, no losses due to area changes as such are included). 

A more realistic inlet model than that of figure 5 would recognize the possibility of 
substantial internal contraction, yielding low Mj, with little loss in pressure recovery. 

In this situation, it is possible that high e's and converging combustors would not appear 
desirable. 


Flight Mach number, Mg 
(b) Burner exit area to inlet area ratio. 
Figure 10. - Supercritical ramjet. 


MIXING 

The preceding section considered the problem of heat addition to a single supersonic 
flow through a nonconstant area duct with applications to the supersonic combustion ram- 
jet. The problem now considered is one of mixing two streams in a nonconstant area 
duct. This applies, for example, to submerging the exhaust of a rocket in a secondary 
flow of air in a ducted rocket arrangement, which results in an increase of specific im- 
pulse due to an exchange of thermal and mechanical energy. A further benefit may be ob- 
tained by a subsequent addition of heat to the mixed flow. 

The algebraic treatment of the three conservation equations, in their integrated form 
for the two flows undergoing complete mixing, leads to complex expressions for the rela- 
tions between certain variables of interest because of the existence of so many parame- 
ters. One of the more tangible relations of mixing two flows is the change of magni- 
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tude of the Mach number of the resultant mixed 
flow M m , which is due to the changes of the 
secondary flow initial Mach number M^ (see 
sketch (f)). Once the character of those vari- 
ations of M m and the bounds on M, for the 
existence of solutions that circumscribe the re- 
gions of interest have been ascertained, it is of 
prime concern to seek to achieve as high a stag- 
nation pressure buildup in the mixed flow as 
possible in those regions of validity. 


Bounds on the Solutions 

With a fixed primary flow, the relation between M m and Mj will depend on the 
mass ratio m^, the stagnation enthalpies ratio A, the static pressure ratio Pj/P^, and 
the power index e when perfect mixing is assumed (see eq. (30) in the appendix). It is 
of interest to first examine how the value of Mj influences the solution for M m . 



< 9 > 



A 

(h) 
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Sketches (g) and (h ) classify the interdependence of M m , M, and A for a specific 
case of c = 1. The three regions of the diagrams correspond to three ranges (low, me- 
dium, and high) of values for the stagnation enthalpy ratio A. In sketch (g), for the low 
and high ranges of A, the shaded regions denote the values of for which there are 
no solutions for M m . Above this shaded region, values of Mj result in an M m that is 
greater than M c , while along the curves separating the two regions the values of 
result in an M that is equal to M . Below the shaded regions, values of M 1 result 
in an M m which is less than M c , and again the bounding curve denotes values of Mj 
that result in an M m equal to M c . For the low range of A, as A increases the ex- 
cluded range of Mj decreases, and finally vanishes at a point M ^ = M m = M c> In 
sketch (h) it is seen that above the shaded region values of M m correspond to values 
of Mj which are greater than M c , while below the shaded region values of M m cor- 
respond to values of which are less than M c> The lines bounding the shaded region 
denote values of M m for which Mj is equal to M c * The excluded region occupies the 
medium range of A. 

The situation can be elucidated further by the familiar N function consideration as 


shown in sketch (i). The relation between N, and N in terms of m 1 and A has to be 
' 1 m 1 

solved to obtain M m< The shaded portions of the 



curve for Nj will produce various values for M m 
until the limiting Mj values are reached, which re- 
sults in M„, at the critical value. A value of M-. be- 
m 1 

tween the limiting points will produce no solution. 
This situation corresponds to the low range of A in 
sketch (g). 

From a return to sketch (h) and an observation 
of the medium magnitude of A, it is noticed that a 



large region of M m values is excluded. This means 
that whatever value for Mj is prescribed, can- 
not assume a large spectrum of values. The corre- 
sponding N function is shown in sketch (j). 

Regardless of the value of Mp M m can be 
found only along the shaded portions of the curve. 

The critical value for Mj results in limiting M m ’s. 

It becomes obvious from the previous discussion 
that the diffusion Mach number Mj delivered by the 
inlet into the mixer is of critical importance, once 
exit conditions have been prescribed, since certain 


tj) 


regions are excluded. Therefore, careful analysis 
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must be done for a fixed configuration flying a changing flight path. 

In order to evaluate overall engine performance, the stagnation pressure buildup 
achieved during mixing and its dependence on Mj must be considered. 


Stagnation Pressure Ratio 


With a fixed primary flow, the stagnation pressure ratio of the mixed flow, besides 
depending in large measure on the secondary flow Mach number Mj, will depend also on 
the mass ratio m^, stagnation enthalpy ratio A, and e (eq. (35) in the appendix). A 
typical graphical representation of the stagnation pressure ratio P^/P° variation 
against Mj for fixed values of the rest of the parameters is shown in sketch (k). 

The existence of two distinct 


Given: 


e - constant 
mj - constant 
A * constant 



Mi 


(k> 




curves is explained by the duality of 

the solution for the mixed flow Mach 

number M . which is obtained 
m’ 

through the evaluation of the N func- 
tion. The situation corresponds to 
the shaded region of sketch (h). 

The two parts of sketch (£) show 
the solutions with their associated 
branches of the curves. These se- 
lected values for M . like I and n 
m’ 

that occur on the same sides of the 
curves as the given M^, are referred 
to as the regular solutions. The two 
on the other sides of the curves will 
be called crossover solutions. In the 
familiar constant- area duct cases 
(e = 1), where only a subsonic Mj is 
allowed, the solution I is elaborated 
and the solution ID, called a super- 
sonic one, is often disregarded. 

A somewhat different stagnation 
pressure ratio diagram (shown in 
sketch (m)) may be arrived at by re- 
collection of the remarks made con- 
cerning the existence of the bounds 
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0 1 2 3 4 5 6 7 

Diffusion Mach number, Mj 


Figure 11. - Regular solutions (I and II) for stagnation pressure ratio at various 
stagnation enthalpy ratios. Mass ratio, 3; power index, 1.0; static pressure 
ratio, Pj/Pj = 1. (See sketch (h) for explanation of gap in curve at A a 0. 10. ) 
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on the solutions for the mixed flow. A spectrum of values of does not give any solu- 
tion for the mixed flow. This corresponds to the situation in sketch (g) (the shaded re- 
gions). In figure 11 the regular solutions (I and D) for the stagnation pressure ratio at 
various A* s have been plotted against for m^ = 3 and e = 1.0. The opposite effects 
of A in subcritical and supercritical mixing are noticed. The gap in the curve for 
A = 0. 10 corresponds to the shaded region of sketch (g). 

It is seen that the selection of is of crucial importance for the performance of 
the mixer, with the value of Mj nearest to unity always yielding the best stagnation pres- 
sure ratio. 

Figure 12 shows the variation of the stagnation pressure ratio of the solutions at con- 
stant M 1 and A against e. The mild increase of the parameter with the increase of e 

is observed. Finally, in figure 13 the dependence of 
the maximum obtainable stagnation pressure ratio 
on the mass ratio is found. As to be expected, the 
smaller the amount of low-energy secondary air, 
the greater the increase in its stagnation pressure. 


Mixing and Burning 


Higher specific impulses may be possible if the 
function of momentum transfer mixing is separated 
from the function of combustion by provision of two 
separate chambers in the duct, one being the mixer 
and the other being the burner. Beside the advantage 
due to heat addition taking place at elevated pressure, 



Power index, £ 

Figure 12. - Variation of stagnation pressure 
ratio with power index. Mass ratio, 3; 
diffusion Mach number, 1.0; stagnation 
enthalpy ratio, A = 0. 77. 



Figure 13. - Dependence of maximum stagnation pressure ratio on mass ratio. Power 
index, 1; stagnation enthalpy ratio, A = 0. 77. 
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another advantage is the opening of scope for the greatest possible stagnation pressure 
buildup in the mixer and the least possible loss of it in the burner, or an optimum com- 
promise for the maximization of the product. It may be expected that the two separate 
chambers, each with its own optimum geometry (in other words, two in number geome- 
tries), show a drastic improvement over the constraint of one single geometry associated 
with a single-chamber concept. 

Ideally, the gross jet thrust per unit of primary mass flow of the engine is given by 

h° — h° — 
h j + q + -L — i— 



whether it is one- or two-chamber variant. 

The difference will be in the value of the stagnation enthalpy of primary jet h?, which 
in the case of the single -chamber case will be quite smaller for limited oxidizer-fuel 
ratios. With the single chamber, the oxidizer-fuel ratio of the primary jet is below the 
stoichiometric ratio, thus yielding a comparatively small h°. It is assumed that the 
fuel-rich primary is the only source of fuel for the secondary flow. 

As to the pressure term, what remains to be mentioned is the comparative difficulty 
of optimizing it in the single- chamber concept. Sketch (n) shows the typical relation be- 
tween M m and plotted against A for the 
single chamber. Large exclusion areas are noted. 
A large spectrum of Mj values produces no solu- 
tion. 

To compare conclusively the virtues of "mix- 
ing followed by burning" to "mixing while burn- 
ing, " an examination of the corresponding stagna- 
tion pressure ratios is in order. 

For every diffusion Mach number Mj and a 
given stagnation enthalpy ratio A, the perform- 
ance of the optimum e mixer duct configuration coupled to an optimum supercritical 
burner (which as seen previously must be a critical e duct) is juxtaposed with the opti- 
mum mixer burner in figure 14. A marked difference is noted. The gap in the lower 
curve is explained by the lack of the solutions for a spectrum of Mp as shown in 
sketch (n). A significant advantage in magnitude of the stagnation pressure and in the 
range of its achievable conditions is seen in the case of separate mixing and burning. 

This advantage prevails in most of the cases. 
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Figure 14. - Stagnation pressure ratios for "mixing followed by burning" and 
"mixing while burning". Mass ratio, 1; power index of mixing duct, 1. O f 
flight Mach number, 10. 


TYPICAL ENGINE PERFORMANCE 

The separate mixer and burner systems need more examination in terms of the over- 
all performance that includes an inlet and a nozzle. 

If for the nozzle a full expansion to the ambient pressure and a given mass ratio m^ 
for the whole engine are assumed, it is necessary to correlate the three stagnation pres- 
sure ratios achieved in induction, mixing, and stoichiometric hydrogen burning at any 
flight Mach number Mq. 

As seen in figure 5 (p. 11) the stagnation pressure ratio of induction is related di- 
rectly to the diffusion Mach number Mj. This variable in turn, as it has been seen in 
figure 11 (p. 18), affects the stagnation pressure ratio of mixing. 
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Stagnation pressure ratio, P^/pJ 



Figure 15. - Stagnation pressure ratios of induction, mixing, and burning. Mass ratio, 
3.0; power index of mixing duct, 1.0; power index of burner, critical; flight Mach n urn- 
ber, 5; separate mixing and burning. 



Figure 16. - Variation of specific impulse with thrust coefficient. Flight 
Mach number, 5; altitude, 70 000 feet. 
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In figure 15, the stagnation pressure ratios of induction, mixing, and burning for the 
case of Mq = 5 and rhj = 3 are plotted against the diffusion Mach number for a 
single e of mixing duct e m = 1. The primary jet is a stoichiometrically burnt hydrogen 
in oxygen, and the static-pressures ratio Pj/Pj is unity. Whatever the resultant mixing 
Mach number M m , the critical burning follows, and its stagnation pressure ratio is also 
shown in figure 15. Repetition of the same procedure for several e m leads to the best 
obtainable overall pressure ratio that can be obtained for the selected case. 

Relative to its specific impulse and its thrust coefficient, the hybrid engine consid- 
ered occupies an intermediate position between pure ramjet and a rocket. The ramjet 
has a high specific impulse but a relatively low thrust-to-weight ratio, and the chemical 
rocket has substantially opposite characteristics. Where exactly a given hybrid engine 
is located in the wide spectrum of both mentioned parameters ought to be implied by the 


entrainment ratio. 

This spectrum, delineated by the ramjet and the rocket at both extremes, is shown in 
figure 16, where the specific impulse is plotted against the thrust coefficient for the flight 



Figure 17. - Variation of specific impulses with flight 
Mach number. 


Mach number Mq = 5 and the altitude of 
70 000 feet. There is an order- of- magnitude 
variation in both propulsion characteristics, 
and the actual choice of the operating point 
will ensue from the integration of the pro- 
pulsion and the vehicle for any particular 
mission. 

In order to have some insight into the 
performance capability of a booster with a 
rocket-ramjet propulsion, the specific im- 
pulse of the engines against the flight Mach 
number has been plotted in figure 17. The 
area enclosed by the curves of the ramjet and 
the rocket maps out the capability of the hy- 
brid engines. One of those engines with an 
entrainment ratio rrij = 5, whose vehicle tra- 
jectory is shown in figure 18, has been plot- 
ted. Its performance has been optimized at 
every flight Mach number by the method 
shown earlier for Mq = 5. Figure 19 shows 
the variation of some important parameters 
of this particular hybrid engine. It is seen 
that a great measure of flexibility would be 
required from the mixer and the burner. 
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CONCLUDING REMARKS 


Based on the assumption of certain idealized conditions, the analysis of nonconstant 
area combustion and mixing has been carried out. The results of this analysis indicate 
that the duct shape is an important consideration in the evaluation of the performance of 
the burner or the mixer. The area ratio as a variable, besides broadening the scope of 
the possible flow solutions, leads to the choice of a suitable and optimum duct for any 
boundary conditions. 

Special attention has to be given to the magnitude of the diffusion Mach number of the 
secondary flow duct entrance. This Mach number affects greatly the performance of the 
duct. 

The nonconstant area combustion in a supersonic combustion ramjet results in a 
better specific impulse. Also, if the separation is feasible, the separate mixing and 
burning, each with its appropriate duct configuration, leads to better performance than 
simultaneous mixing and burning in the case of the augmented rocket. 

Admittedly, the varying flight conditions could in general demand a varying configura- 
tion or many different fixed configurations. One of those configurations, however, will do 
better for the whole flight spectrum than the rest of them; thus, a fixed- geo me try engine 
design is not precluded. The best geometry can be ascertained by the derived methods of 
this report. In general, the insight into nonconstant area processes afforded by this re- 
port should be helpful in the design of actual engines. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 3, 1966, 

126-15-09-08-22. 


25 



APPENDIX - EQUATIONS OF BURNING AND MIXING 


Various relations are derived, starting with the standard one-dimensional treatment 
of the three conservation equations as modified by Crocco's device for the use of single- 
flow heating. 


Heat Addition to a Single Flow 



The conservation of mass equation is (see sketch (o)) 

mj = m 2 (1) 


or since 


m = P Av = ^\€yir» PAMV ^ ° PAyFv ^ 
RT ^ 


where 


it will assume the form 


F = M 



y ~ 

2 


1 


M 


P 1 A 1 F 1 P 2 A 2 F 2 



( 2 ) 


if r 1 =r 2 

From the conservation of momentum. 


1 

e 


( m 2 v 2 - m l v p ~ P 1 A 1 " P 2 A 2 



If it is assumed that the Crocco relation holds. 


( 3 ) 
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P (. A V/( 1-e ) 


P 1 \ A 1 


P dA 


■M 


e/(l-e) 


dA = (1 - e) (P 2 A 2 - PjAp 


Substituting this into equation (3) gives 


m 9 V 9 m 1 V 1 

— — + eP 9 A 9 =-i-i + eP t A t 
g 2 2 g 11 


Since 


V = MVgyRT 


m _ PAMy\/g 


Vy RT 


the equation can be transformed into 


P 2 A 2 (€ + r 2 M 2) = PjAj ( £ + yjM^) 


P 2 A 2 G 2 “ P 1 A 1 G 1 


where 


G = 6 + yM 


The conservation of energy equation is 
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( 5 ) 


h j + q = hg 


In equation (2) F can be replaced by GN where 




M 


e + yM 


2 


( 6 ) 


Then 


P 1 A 1 G 1 N 1 P 2 A 2 G 2 N 2 



which on account of equation (4) becomes 





Assuming h^Aj = Tg/Tj and substituting equation (5) into it yield 



This equation is used to obtain M 2 , the postcombustion Mach number. The critical value 
of Mg is derived from dN/dM = 0, which from equation (6) becomes 
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Stagnation Pressure Ratio of Burning 


From equation (4), 



A 1 G 1 

A 2 G 2 


which on the substitution of Crocco's relation 



/ A \€/(i-e) 
I — I 
\ A V 


becomes 


Also, 


Then, 


where 


or 



( 9 ) 



( 10 ) 



X = 1 + 


-Am 2 
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where 


G* = °! 


(y- i)/r 


x 


Specific Impulse 

Velocity of the jet at point 3 after expansion follows from sketch (p) and equation (5): 



s 

(p) 


P? p2 p? 


X 


y/(y-i) _ 3 _ *3 _ v-y/ (y _ 1) _1 _2 __3 


"o 


p° p° p° 
P 0 F 1 P 2 


it follows that 



<.) 


X ° \ P ? P 2 P ?/ 


( 12 ) 


and then 
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( 13 ) 


V 3- V 0 


Critical Epsilon 

Substituting equation (8) into (6) gives the value of N function at the critical point: 


N 2 = — I 

4ey j _ e y - 1 


Deriving an explicit expression for M from equation (6) yields 


1 - 2eyN 2 ±l/ 1 - 4N 2 ey fl - ey 


2y 2 /n 2 - 


When the numerator and denominator are multiplied by 


1 - 2eyN 2 Tl/l - 4N 2 ey[l - ey 21^-1 


it follows that 


4e 2 y 2 N 2 (N 2 - 

\ 2y 2 , 


2y 2 ^N 2 - 1 " 2eyN 2 *y 1 - 4N 2 ey^l - ey 


Using this equation and equation (14) gives 
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Multiplying both sides of the equation by N 2 yields 

V 



However, 


N 2 m2xG c 

N 2 M 2 X G 2 
c c c 


which means that 


N 2 M c G 2 . N 2 M c M 4 N c 


N 2 M 2 G 2 
c c 


N 2 M 2 M 4 N 4 
c c 



n 2 m 2 


Hence, it follows that 
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2 


or 


where 


Hence, 




sin a 


M 2 _ /m 2 \ sin a 
F y F y 1 ± cos a 


(15) 


Substituting Mj on the left hand side of equation (15) and 


N, 


sin a = 


N„ 



on the right hand side will give the value of M . Then from equation (8) the critical e 

V 

is obtained. 

Hence this equation makes it possible to obtain directly the duct configuration that 
will secure a critical burning for a given Mp heat condition, and stagnation enthalpy. 
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Mixing Equations 

An algebraic treatment of the three conservation equations is now presented (see 
sketch (f), p. 15). 

First, from the conservation of mass, 

m l + m r m m 

In a manner similar to the way equation (2) was obtained, this can be put in the form 

P l A l y l F l P )Vj F i _ P m A m>-m F m ( 

VS? 

Second, from the conservation of momentum, 


i • r m 

- f m m v m - m l v l - m j v j] = P 1 A 1 ♦ Pj A j - P m A m t J P dA 


(17) 


If it is assumed that the Crocco relation holds, 


ji = ( a V ^ 1-e * 

P 1 \ A 1 + A j/ 


then 


/ m /»m / \e/(l e ) 

PdA= / Fl y^y dA=( 1 - £)[PmAm ' p i (A i tA i )I (I8) 


Substituting equation (18) into equation (17) gives 


l t m m v m - m l v l - ' P 1 A 1 * PjAj - PjAq - + sPjfA, + Aj ) 
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m_V m.V, 


m.V. 


eP A + = i + eP 1 A 1 + — U + ePjA, + P,A. - P t A, 

mm _ _ 11 g 1 J 3 ] 1 J 


g g 


n^V* m.V. 

+ eP.A- + — L-I + ePjA. + P^, - PjA. + eP^ - eP.A. 

g g J JJ J JJ J J 


m.V. m.V 

+ eP^j + —LJ. + ePjA. + (P. - P t )A.(l - e) 

g J.A g JJ J A J 


m l v i m i v i 

i + €Pj A- + — L-i + P. A. 

g g J J 



€ + [1 - -±\ (1 - e) 


m,V 


m.V. 


* * + ePjAj + — Li +e p a. 

A g J J J 


where 


= c + (l - — 1(1 - e) 


Since V = M VgyRT and m = PAMy /VyRT, the previous equation can be put in the 
following form: 


P m A m( £ + >-m M D ' P l A l( £ + y l M !) * P Mi * y i M f) 


or 


P m A m G m = P 1 A 1 G 1 + P j A i G J 


( 19 ) 


where 


G = e + yM 2 

From this form of the conservation of momentum equation, 
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but 



A i G i + 7' A ) G j 


m m 


A m = < A 1 + A j> A 


m 


where 



m 


a i + A j 


Hence, 



Pi 

A 1 G, + A.G. 

1 1 p ] 3 G 


(A* + A.)A* G 
v 1 j' m m 


av 

A* G 
mm 


where 


P. _ P* 

A, G 1 + — A^G- A 1 G 1 + — G. 
1 1 Pi 3 3 1 1 Pj J 

Aj + Aj Aj + 1 


From Crocco's relation, 


so 



e/(l-e) 




ve/(l-e) 


V p i 


( 20 ) 


( 21 ) 
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Hence 


and 



Third, since 

m _ PAyFVg^ _ PAyF\^ 

Jyl 2^1 h° 

f y 

and equation (1 9) can be put in the form 


F, F. F 

p ‘ A i s: + p i*j i ' P ™ A ™ if 

1 ] m 


one can substitute equation (24) into this equation with the result 



From conservation of energy, 

m, h? + m.h? = mh° 
11 j 3 mm 


and this equation can be put into the form 


( 22 ) 


(23) 


(24) 


(25) 


(26) 
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( 27 ) 


1 -f m 1 ^ — 

h° 1 h? 1 + m. A 

m j_ _ 1 

h? 1 + rri| 1 + nij 

Substituting equation (27) into equation (25) gives 



This is the fundamental equation for evaluation of mixing duct end Mach number M m 
when the entrainment ratio m^, stagnation enthalpy ratio of two flows A, and Mach num- 
bers of both flows are given. By simple transformation, equation (28) can be put in the 
form (with the assumption y j = y m = y^) 


or 


where 



A. 

N m 


m 




^/(l + m,) (1 + m, A) + m,) (1 + m, A) 


A 

N* 


(29) 


(30) 


- T P 1 F 1 1 
= A 1 

1 p i F iv£ 
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Stagnation Pressure Ratio During Mixing 


The stagnation pressure ratio of the mixing process is now evaluated. 
From conservation of momentum equation (19) 


P mV F m _ P 1 A 1 F 1 + P j A j F j 


it follows that 



N m 


N 1 

N J 

P m 

N m 

/ A 1 F 1 

+ Zl 

Vi) 

P 1 

**i 

S 

U 

p i 

"J 


N 

m 

/ A 1 F 1 

+ H 

A. F. 

J A 1 


A m F m 

V N i 

P 1 



A 1 F 1 


Solving equation (24) for PAF and substituting in the previous equation yield 

\ 

1 


P m _ N m A l F l 
P 1 A m F m 


1 - + — L 


\N 


1 m 




if y - constant 


A 1 F 1 / N m 


N_ 


m 


^m F m V N 1 ' T ™ 


Njirij y/Kj 


(31) 


Now since 


4. A i A i tA i 


Am A 1 + A j A m 


1 + A, A 


-f n 

1 m 


(32) 


substituting equation (32) into equation (31) in turn results in 


p 

* m 


A 1 F 1 / N m 


N_ 


m 


m 


P 1 1 + A t F m 


N« XT — 


Njm 1 VaJ 


(33) 
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Crocco's relation gives 


-£ = (A* ) e/(1 “ e) 
p m' 


From this 


and 







Substituting equation (34) into equation (33) results in 



Then the stagnation pressure ratio will be given by 



Hence, since N/F = 1/G, 
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(y-i )/y 


e(y-l)/y 


m 


a + A i 


i+- 


Nl \ e(y " 1)/y °i 


N^m^VA^ 


g: 


m 


(35) 


where, as in equation (11), 


* _ G 

G — 


Burning After Mixing 


Conservation of energy is 


m m h m + “l* = m m h 2 


where 


m = m. + m. 
m 1 j 


Hence, 


(mj + l)h^ + mjq = (n^ + ljhg 


(36) 


During mixing, the conservation equation was 


h? + nijhj = (m^ + l)hj 


o 

m 


(37) 


Substituting equation (37) into equation (36) gives 


h° + rajhj + = (nij + l)h° 


(38) 


Then 
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2 


_ T 2 _ + m l h l + m l q _ 


„ 1 + m 1 A + m 1 -S- — 

1 1 h% 1 + nij A 


IN, 




m 


h? + m^hj 


1 + mj A 


1 + mjA 


(39) 



where 


oo a 

A = A + -^- 


The velocity that can be developed at point 3 (see 
sketch (q)) will be given by 


2gJ 


= h 3- h 3 = h 3 



i,° 



o — 


hj + mjhj + m^q 
1 +IU, 



Then gross jet thrust per unit primary mass flow is 


(mj + 1)V 3 = V2gJ (1 + m 


/h? + m h° + m q 

JV - — 

f 1 + m^ 



and 
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Burning While Mixing Process 


The primary jet here is stoichiometric and the fuel has to be supplied for the stoichi- 
metric combustion of the entrained secondary flow. The mass of fuel is considered negli- 
gible as compared to m^ + im. 

The conservation of energy equation can be stated as 


or 


mjhj + m^h? + m^q = nighg 
m^hj + h? + m^q = (m^ + l)h£ 


This is the same as equation (38); hence, the same equation for the jet thrust as shown 
in (40) is obtained. The specific impulse is 

(mj + 1)V 3 - mjVQ 


The primary jet here is fuel- rich and constitutes the only source of fuel for the secondary 
flow. 
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